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1 Introduction 



Among the models put forward for creation of baryon asymmetry, those based on elec- 
troweak baryon number violation during the electroweak phase transition (EWPT) 
stand out because the physics involved in the relevant processes is either being tested, or 
can soon be tested in accelerator experiments. A good example is the sphaleron wash-out 
constraint: to maintain an asymmetry created in EWPT, the baryon number violation 
must turn off in the broken phase. Whether this is the case, depends on the strength 
of the transition f^, which again depends on the model parameters, and can be studied 
using either effective potential or lattice methods. In particular it has been shown |^ 
that, given the current bound on the Higgs mass, the transition in the Minimal Stan- 
dard Model (MSM) does not satisfy this constraint. So, to explain the origin of baryons 
one is led to consider extensions of the MSM, and perhaps the most natural is to in- 
troduce supersymmetry (SUSY). Among SUSY extensions of MSM, already the simplest 
Minimal Supersymmetric Standard Model (MSSM) and its next to minimal extension 
(NMSSM) 1^ contain additional scalars which can strengthen the phase transition as 
required for baryogenesis. 

Because of the out-of-equilibrium requirement 0, baryons can only be created within 
the moving phase transition walls or in their immediate vicinity during EWPT. It is no 
surprise then that the ratio of the wall width £„aii to the other characteristic scales, such 
as the scattering lengths 4catt ~ 1/r and de Broglie wave length = 1/p ~ '^/T of 
the particles involved, strongly influences the mechanism of baryogenesis. The wall width 
also controls the CP-violation, which in general is introduced through complex phases 
in the mass matrices. When mass matrix elements depend on the higgs condensate, the 
phases of mass eigenstates can change over the wall profile influencing the propagation 
of particles and antiparticles in different ways. If the wall is wide in comparison with 
the typical de Broghe wave length f^aii 4b, these effects should be computable using 
semiclassical (SC) methods 0]. 

In supersymmetric models the bubble walls are found to be quite slow and thick ^, 
so that the semiclassical approach [0] should be applicable. It has been used to compute 
the baryogenesis from chargino sector of MSSM [jlO|, [IT], |12| and in the NMSSM [|l3l. How- 
ever, the semiclassical approach is based on the assumption that the plasma in a spatially 
varying background can be described by a collection of single particle WKB-excitations. 
Despite the obvious intuitive appeal, this assumption was lacking a fundamental proof 
until recently |T^. Proving the validity of the SC-approach was important, in particular 
because many other formulations of the transport equations in the limit of slowly varying 
walls had been but forward, which gave both quantitatively and qualitatively different 
results from the SC predictions Moreover, there was some initial confusion in identi- 



fying the correct physical excitations and their distribution functions even within the SC 
picture J^, |10[, and a fully consistent SC-computation was introduced only in ref. [|lT . 



Here we review our derivation of coUisionless SC-transport equations for fermions with 
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spatially varying CP-violating mass terms |]14| starting from Schwinger-Keldysh equations 
for the dynamical two point Green function (Wiglitman function) C^lujv). Moreover, 
the analysis of [|l4l is extended to include the gauge fields, which give rise to the usual 



Lorenz-force on gauge-charged particles. We also generalize the WKB- analysis of ref. [11 



to include the gauge fields, and show how the gauge degree of freedom is connected to the 



"reparametrization dependence" of canonical variables in the WKB-picture. As in 
we concentrate on the 1+1-dimensional model; the full 3+1 dimensional case will be 
considered in ref. [ |T^ . 

The most remarkable result of our analysis is that to first order in gradients (or 
equivalently in h), the equations of motion admit a spectral decomposition solution for 

in terms of on-shell quasiparticle excitations; this qualitatively proves the validity of 
the above mentioned assumption of single particle picture underlying the SC-approach. 
Moreover, the on-shell momenta set by a dispersion relation derived from the equations of 
motion quantitatively agree with the SC results of |ll| (and also with an earlier Schwinger- 
Dyson equations based calculation [|T7|). Finally, the on-shell distribution functions are 
defined for the quasiparticle excitations, and are shown to obey the usual kinetic (Vlasov) 
equations, with the SC group velocity Vs± = kz/ujs± and the force Fs± = ujs±dvs±/dt. 

However, the equations of motion only admit a spectral decomposition solution to first 
order in h. At higher orders more complicated structures arise and the SC-picture of a 
plasma consisting of single particle excitations breaks down. Quite fortuitiously for the 
EWBG-program this is all we need because, as we shall see, the dominant CP-violation 
terms appear precisely at order h. 

This contribution is organized as follows. In section 2 we derive the SC WKB-equations 
for a single Dirac fermion with a spatially varying complex mass term including the gauge 
field. In section 3 we derive the SC kinetic equation from Schwinger-Keldysh formalism, 
and in section 5 we study vector and axial vector divergences, and show how one sets up 
moment equations for practical baryogenesis computations. Finally, section 5 contains a 
discussion and summary. 



2 WKB method 

We will here derive the semiclassical dispersion relation for a fermionic field with a complex 
spatially varying mass term. More precisely, we take our system to be described by the 
effective lagrangian 

£ = iij{p+ie4)i^ - i^imipR - ipRm*i)L + Ant , (1) 
where is the [/(l)-gauge field, £int contains (here unspecified) interactions and 

■m{z) = niR^z) +imj{z) = \m{z)\e'^^^^ (2) 
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is a mass term arising from an interaction with some CP- violating scalar field condensate 
with a planar symmetry. Planar symmetry effectively reduces the problem to 1+1 dimen- 
sionsQ (by use of boosts), and to keep the notation simple and consistent with the rest of 
the paper, we will assume this to be strictly the case here. 

The basic assumption in the semiclassical treatment is that the plasma can be approx- 
imated by a collection of single particle excitations, whose dispersive properties can be 
studied by employing the WKB-methods. The lagrangian (|lD implies the following Dirac 
equation for the spinor wave function %l){z,t): 

(jp — e ^—tur — i'y^mj^ ip = 0. (3) 

We are interested in stationary solutions of the form ip{t, z) = e~*^''V(^o; z)- Eq- ® then 
becomes 

(7°fi:o + i-f^{d, - ieA,) - thr - i-f^mj) iP{ko; z) = 0, (4) 

where kq = ko — cAq. Here spin is also a good quantum number and we can write the 
wave function as a direct product of chiral and spin degrees of freedom: 

^^[rI)^x', (5) 

where are the eigenspinors of the spin operator a^. Eq. @) then splits to two coupled 
equations or the left and right chiral wave functions: 

(kq — esAz — isdz)Ls = mRg 

[kq + esAz + isdz)Rs = m*Ls, (6) 

One can use, say, the first of equations (H) to eliminate Rs from the second one, to 
obtain a second order equation for Lg. That equation can be analyzed by making the 
WKB-parametrization for L^: 

Ls=Use'r^''^^' . (7) 
In this way one finds the following two real valued equations 



u'' \m\' u' 



ni- K\-\mY + {skq + K^)e' + — - = 



u' , , Iml' 



4 + (2fi:3-^')— -(s'^o + '^s)^ = 0, (9) 



where = ks — eA^. These equations are exactly analogous to the ones obtained in |TT 



with the sole exception that the bare momentum of |]rT|, has here been replaced by the 
one including the gauge field: 

kfj, ^ i^fjL = kfj_ — eAfj^. (10) 



In partiuclar the gauge field components A^^y can be taken to vanish due to planar symmetry. 
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In absence of the gauge field, the WKB-equations of ||TT[ were variant under phase re- 
definitions of the ip-fie\d, which led to the "reparametrization variance" of the canonical 
momentum. Equations (H) and (P) are clearly invariant under complete gauge transfor- 



mations ^p — > e^^'ijj and — {i/e)d^a, and the reparametrization variance of [11] 

is here seen to correspond to the usual freedom of adding any four divergence of a scalar 
function 9^0 to the gauge potential A^. 

While complex in appearance, equations (|[]^) have the advantage of being readily 
solved iteratively in gradient expansion. This part of the analysis is analogous to one 



shown in |11], and hence we only quote the final results. To the first order in gradients 



the dispersion relation reads 

{sko+po)e' 



Ks = h- eAf, = po + — , (11) 



where po = sign{ps)ykQ — |mp and s^= 1(— 1) for particles (antiparticles). Equation 
(pup can be inverted to give the energy in terms of the momentum variable ^3: 



ft:o =^ y (/t3 - s^—f + |m|2 - s^. (12) 

Now it is easy to compute the physical momentum, which is defined as usual, in terms of 
the group velocity kz = i^oVg = i^o{dK^Ko)z- After some algebra one finds 

s\m\'^9' 

^^=^° + >-2^- ^^'^ 
Finally, the physical force is defined as = k^. The only difference between the compu- 



tation here and in ref. |]TT| is that in the presence of an electric field the energy is not 
conserved; indeed after some algebra one finds the expected result = eVgE^, where 
is the electric field in the z-direction. Similarly, the classical force is found to be 

Imp' s(\m\^9')' ( s\m\^d'\ , , 

Fz = kovg + novg = + s^-^^ + eEz (^1 - s^-^ j . (14) 

We have here computed all expressions to the first order in gauge field and to the second 
order in gradients. If the electric field Ez is a self-consistent field created by the charge 
separation due to CP- violating force in (|T^, it is formally a second order term and the 
gradient correction multiplying eEz should be dropped. According to the semiclassical as- 
sumption the Vlasov equation for the phase space density of coUisionless WKB-excitations 
can now be written as 

dtfs + Vgdzfs + Fzd,J, = 0. (15) 

The above computation is improvement upon ref. |TT| in that we include also the self 
consistent electric field which is necessarily generated by the charge separation induced 
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by the wall. It also better explains the phase reparametrization variance observed in 
connection with the canoncial momentum. The WKB-derivation however depends on the 
important assumption that the plasma in a spatially varying background can be described 
in a single particle WKB-picture. To determine whether this really is so, and to which 
accuracy, one has to go beyond the semiclassical approximation. This is what we do in 
the next section. 



3 Schwinger-Keldysh method 

The statistical properties of an out-of-equilibrium system can be described by the two 
point Wightman function 

G^piu,v) = z{Uu)e^^Syy^Mv)). (16) 

where (■) denotes the expectation value with respect to the initial state and the path for 
the gauge flux is defined as 



ru rl/2 

/ dy-A= dsr^'AJx + sr), (17) 

Jv J-1/2 



where r = u — v denotes the relative and x = {u + v)/2 the center-of-mass coordinate. 
The definition (|16D is a gauge-invariant generalization |19[ of the 2-point function used 
14|. corresponds to the off-diagonal part of the fermionic two-point function in 



m 



the Schwinger-Keldysh formalism |^0|], which indeed is the method of choice to derive the 
equations of motion for including interactions. However, for noninteracting system 
studied here, all one needs is the familiar Dirac equation. Using Eq. (0) one finds: 
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miuYPL - m(u)PR + e-i^—- / dx ■ A] G<{u,v) = 0. 

ou^ Jv I 



In order to study equation (|T8|) in gradient expansion we performa a Wigner transform 
of to the mixed representation: 



G<(x, k) = I d^r e'''-'G<{x + r/2, x - r/2), (19) 



where r and x are defined as in (p!7|). In the mixed representation the equation of motion 
becomes 

(j^ik, x) + '-I^{k, x) - mnix - '-dk) - imiix - G< {k, x) = 0. (20) 



where 



rriRjix - -dk) = mR,7(a;)e-^^--^'= (21) 
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and the gauge generalized derivative and momentum terms are (in agreement with Zhuang 
and Heinz 

/•+1/2 

J — 1/2 
/•+1/2 

n^(x, /c) = k^-ie dssF^^{x - isdk)dk,. (22) 

J-l/2 

Above equations of motion are manifestly gauge invariant, and one can in fact always 
rewrite the F^i,-terms in any of the above formulas in terms of physical electric and 
magnetic fields. Using = Ei and Fij = —eij^B^. one finds 

YF^^du^ = -^E .d,^ + ^-Edk,+i-{Bx d,^). (23) 



The crucial advantage of the representation ( |19|) is that it separates the internal fluc- 
tuation scales, described by momenta k, from the external ones which show up as a 
dependence of on the center-of-mass coordinate x, and thus gives us the chance of 
exploiting possible hierarchies between these scales. Indeed, when the external fields vary 
slowly, one can truncate the Taylor series operator expansions of the mass and gauge field 
operators appearing in (pi|-PB|) to some finite order in gradients (or equivalently in powers 
of h). To the lowest order we get 

Di,{x,k) ~ - e Fi,^{x)dk, 

U^{x,k) ~ k^. (24) 

In the planar symmetric case of an electroweak bubble wall, there can be no magnetic 
field in the plasma frame. Also can only depend on the spatial coordinate orthogonal 
to the wall, z = x^. Continuing to work in the strictly 1+1-dimensional case like in the 
previous section, we also drop the /c||-term. With these simplifications the equation (|20|) 
reduces to 

{ko + kzl^l^ - mo7° + ^^57%^) i7°G< = 0, (25) 
where we used the following shorthand notations 

ko = ko + -{dt + eEzdkJ 

k = kz - '-{dz + eEzdko) (26) 
where is the electric field in the ^-direction, and 

^0(5) ^ mn(i){z) - -m'^^j^dk,. (27) 
Apart from the electric field Ez appearing in the definitions (EUp , equation (P3|) is identical 



to the one we studied in [14]. Hence the further analysis proceeds in very a similar fashion. 
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we can skip many of the details here. Due to conservation of spin in z-direction can 
be decomposed as follows: 

-il'G< = \{l + sa')®gf, (28) 

where is the usual Pauli matrix. Using this structure, and making the appropriate 
identifications 7° p^, — Z7°7^ — * and —7^ p^, the original 4x4 problem reduces 
to a two-dimensional one: 

{ko - skp^ - P^rfiQ - p^rh^) gf = 0, (29) 

where are the Pauli matrices, which span the remaining chiral indices. The hermitean 
matrices gf can be expanded as 

9f = I {90 + 9tP') , (30) 

where the signs and normalizations are chosen such that gQ measures the number density 
of particles with spin s in phase space. 

Equation (^) still constitutes a set of four coupled complex differential equations for 
components of g^ (a = 0,^). They can be obtaned by multiplying ( p9D successively by 1 
and p* and taking the trace: 



^ofl-o - sh9'3 - ^o9t - rh^g'^ = (31) 

kogi, - sk^gQ - irhog?. + irn^g{ = (32) 

kogl + isk^g^ - mofi-o - im^gl, = (33) 

kog2 - ishgt + im^gl - m^g^ = 0. (34) 



The real parts of these equations contain no time derivatives and hence are called con- 
straint equations (CE). Their very important role is to constrain the available space of 
solutions for the four dynamical, kinetic equations (KE), which correspond to the imagi- 
nary parts of ( pTj - p^ ). 

3.1 Constraint equations 

We first use three out of the four constraint equations to iteratively solve gf in terms of gQ 
and its derivatives. These expressions are then used to eliminate g^ from the last remaining 
CE. The most remarkable property of this equation is that all terms proportional to dk^g^, 
also those involving the field Ez, cancel to first order in gradients (or h). As a result g^ 
satisfies an algebraic equation 

^kl-kl-\m\' + ^Jm\'9'^g^Q = 0, (35) 
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which admits the spectral solution 



^0 = IlT7^^s^(fcoTc^s±). (36) 

Here ns{kQ, kz, z) are some nonsingular functions, the on-shell energies ujs± are specified 
by the roots of equation 

= kl - kl - |m|2 + f |m|2^' = 0, (37) 

and the normalization factor is defined as Zs± = 2;^\dko^1\ko=±u}s±- 

The solution (|3^) is in many ways the most important result of our analysis, because 
it precisely proves the WKB-assumption that the plasma, to a certain accuracy, can be 
described by a collection of single particle on-shell excitations. Below we shall further 
show that also the WKB-dispersion relation and the semiclassical equations of motion of 
ref. |Tl| do arise in this picture. It is important to note however, that the confinement 
to sharp energy shells is only valid up to first order in gradients; beyond first order the 
constraint equation for Qq is no longer algebraic and hence does not admit a spectral 
solution, implying breakdown of the SC picture. 

We end this subsection noting that, when treated to the first order in gradients, 
Eqn. ([37| ) yields the solution 



ujs±=uJoTs 2 , uJo = Jk^^ + \m\'^. (38) 

ZLOq 

Inverting this to give kz in terms of energy, one finds it to be identical with the WKB- 
relation (^) for the physical momentum. 



3.2 Kinetic equations 

We now turn our attention to the kinetic equations. One could choose any one out of 
the four equations to describe the dynamics along with the constraint equations. The 
quantity of most direct interest is Qq however, since it carries information on the particle 
density in phase space. From (pT]) one gets 



{dt + eEzdk,)g'o + s{dz + eEzdko)gl - m'^dk^g{ - m'jdk^g'^ = 0. (39) 
Using constraint equations to eliminate g^ as in the previous section one finds 

' i\mW-eEz(l-'-^]]dk^g^, = 0. (40) 



2/lq \ 
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We have so far used three out of four constraint equations. The fourth constraint (|36|) 
forces the solutions to the single particle on-shell, and is easily accounted for by inserting 
the spectral solution (|36D into (|40|) and then separately integrating over the positive and 
negative frequencies ko. The total divergence term with respect to ko drops and get the 
familiar form for the Vlasov equation 

dtfs± + Vs±dJs± + Fs±dkJs± = 0, (41) 

where fs+ = ns{ujs+, kz, z) and fs- = l — ns{—ujs-, —kz,z) are the particle and antiparticle 
distribution functions with spin s, respectively, the quasiparticle group velocity is given 
by Vs± = kz/ujs± and the spin-dependent and CP- violating semidassical force reads 



in complete analogy with the semidassical result. The physcial significance of ( |42| ) is that 
due to derivative (quantum) corrections the spin degeneracy is lifted at first order in gradi- 
ents and hence particles (antiparticles) of different spin experience different accelerations 
in a spatially varying background. 



4 Continuity equations and sources for BG 

We now investigate how the CP-violating sources appear in the equations for the diver- 
gences of the vector and axial vector currents. First note that any current operator can 
be computed as a trace in the Wigner representation, weighted by G^: 

(V5(x)X^^(x)) = / (0^Tr(-^G<X^). (43) 

For vector current = and for axial one = 7^7^. Using the above expressions 
and constraint equations from (|3T1 - |3^ one can verify that the divergence of the vector 
current is conserved and the axial one is broken by the complex mass as expected: |^ 

d^f = (44) 
d^j^ = 2iMn{^^''i;)-2Mj{i,i;), (45) 

TNote that the sign of the m/j-term in the axial current divergence was incorrect in our original 
work |l4| . That error was compensated however, by a another trivial sign error in the definition of the 
pseudoscalar density in terms of . 
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where (ips'ips) = Jk9i denotes the scalar, and (ipsl^i's) = —i Jk92 ^^e pseudoscalar density, 
and Jk = J cPk/{27iY. If we define the fiuid density^ ns± and fiuid velocity moments as^ 

(2^ 

1 r d^k (K^^ 



= — / 7^ rMo' (46) 

ns+ J+ (27r)^ \ fco/ 



then Eqs. ( ^^^31 ) can be written as: 



dtns± + d^{ns±{vs±)) = 0, (47) 
dt{ns±{vs±)) + djns±{v'',^)) = eE,\m\''lss+ + Ss±, (48) 



where the quantity Ss± is related to the baryogenesis source]^ |[n]|, and is given by the 
average over the semiclassical force (E2D divided by ko: 



Ss± = ~\\m\^'l2s± + \s^s{\m\^e'yi,,±. (49) 

with 

In anticipation of the fact that is a self-consistent field, we dropped the gradient 
corrections to the Lorenz force above. The |mp-suppression in the £^2-term is particular 



to 1 + 1 dimensions, and to our use of velocity moments in fiuid equations as in ref. ||1 1|| . 
In the actual baryogenesis computation one expects that the role of the electric field is 
expected to be taken by the hyperelectric field corresponding to the hypercharge U{1)b 
symmetry. 

From expressions ( ^Tj) and ( ^Hf ) it is obvious that the vector current divergence cor- 
responds to the zeroth and the axial current divergence to the first velocity momentum 
of the kinetic equation (03). The former provides just the usual continuity equation, and 



the nontrivial CP-violating source and the electric field only appear in the latter. This 
connection with the axial current nicely explains why the semiclassical source appeares 
only at first order in moment expansion in earlier treatments 0, 11 . 



The fluid density ns± should not be confused with the phase space density in (^6|). 
^ Here we defined integration over the positive frequencies. To obtain equations for antiparticles one 
integrate over negative frequencies taking into account the appropriate thermal equilibrium limit as in 
derivation of (^). 

^Note that the equations ( ^7[]4^ ) are somewhat formal, and require more work before useful for practi- 
cal computations. In particular the velocity-dependence of the actual source only comes about after one 
expands various quantities in the wall velocity around the equilibrium distribution. For details see ||ll| 
and 
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Equations ( ^^ - |48|) do not close, because the field is not yet determined and since 
the two equations contain three other kinetic unknowns n^-t, {Vs±) and The former 

problem can be handled ||22| by adding the Gauss law 

d^E^ = e'^{ns+ - ris-), (51) 

s 

which completely defines the self-consistent electric field. The latter problem is a standard 
one for the moment expansions and requires a truncation approximation. For example one 
may set — > {vs±)'^ + <7s± ^^e the equilibrium distributions for /s±'s to evaluating 
the variance a^^ = — {vs±.)'^ and the integrals Xps± appearing in the source term 



The first source-term in (|49|) is primarily relevant for phase transition dynamics |g, 
but it also contains a spin dependent and CP-violating part coming from the wave 



function renormalization |]T4|. That term, together with the second source provide the 
charge separation necessary for baryogenesis. To promote equations (^7|-^) and ( pT| ) 
into transport equations for baryogenesis calculations we still need to generalize them to 
include collisions. For a discussion of this problem see the talk of Tomislav Prokopec in 
these proceedings [p^ . 



Let us finally discuss the role of the gauge anomaly in the fluid equations. In the 
1+1-dimensional case the U(l)-theory is anomalous, and in general the anomaly [p^ 



(^/iJs ) anomalous 2,71^^'^^ (^^) 

should be added to the axial vector divergence (^). Because the axial current divergence, 
apart from the anomaly, has been here shown to coincide with the first velocity moment of 



the fluid equations |]14[, we can draw the following corollary: the first two moments of the 
fermionic transport equations, with the anomalous term added to the latter, together with 
the appropriate dynamical Langevin equations of motion for the gauge fields, provides a 
complete and self-consistent method for treating the CP-, P- and fermion number vio- 
lating (baryon number violating) out-of-equilibrium interactions required for electroweak 
baryogenesis. 



5 Discussion and summary 



The question of a first principle derivation of CP-violating fluxes in transport equations 
has been the main theoretical challenge of recent work on electroweak baryogenesis. In 



this talk we reviewed our derivation of the kinetic equations appropriate for EWBG 
in a systematic and controlled gradient expansion from Schwinger-Keldysh equations for 
two point Wightmann function, and extended it to include the gauge fields. We also 
generalized the derivation of WKB-dispersion relations |TT 
of freedom. 



to include the gauge degree 
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We have shown that to first order in h the collisionless kinetic equations for fermions 
allow a single particle description with a distribution function which obeys a familiar 
Vlasov equation where the group velocity and the semiclassical force terms contain all 
relevant quantum information, and in particular the CP-violating terms which source 
baryogenesis. These results agree with the SC equations derived in the WKB-picture in 
Ref. originally developed for EWBG problem in 0]. A novel new feature discovered 
in our approach is that the semiclassical limit is only valid to first order in h. Fortunately 
this is not a serious drawback to the EWBG-program, as it is just the first order term 
that carries the dominant source for baryogenesis. 

For simplicity we have considered only the collisionless limit in 1+1 dimensions. One 
can show [jl6| that generalization to the 3+1 dimensional case does not affect results 
qualitatively, and some aspects of collision terms are considered by Tomislav Prokopec in 
these proceedings. For the lack of space within this presentation we also restricted us to 



the case of a single Dirac fermion. For the case of mixing fermionic fields see ref. [14 
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